Effects of Ti Scattering on Fourier-Transformed Inelastic Tunneling Spectra in High-Tc 

Cuprates with Bosonic Modes 
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We study the n-impurity induced q-space pattern of the energy derivative local density of states 
(LDOS) in a d-wave superconductor. We are motivated in part by the recent scanning tunneling 
microscopy (STM) observation of strong gap inhomogeneity with weak charge density variation 
in Bi2Sr2CaCu208+6 (BSCCO). The hypothesis is that the gap inhomogeneity might be triggered 
by the disorder in pair potential. We focus on the effects of electron coupling to various bosonic 
modes, at the mode energy shifted by the d-wave superconducting gap. The pattern due to a highly 
anisotropic coupling of electrons to the Big phonon mode is similar to preliminary results from 
the Fourier transformed inelastic electron tunneling spectroscopy (FT-IETS) STM experiment in 
BSCCO. We discuss the implications of our results in the context of band renormalization effects 
seen in the ARPES experiments, and suggest means to further explore the electron-boson coupling 
in the high- Tc cuprates. 

PACS numbers: 74.25.Jb, 74.50.+r, 74.20.-z, 73.20.Hb 



The extent to which collective excitations of high- Tc 
cuprates are manifested in their single particle spectra is 
a long standing issue. The band renormalization effects, 
seen in the ARPES Q (as well as in the break junction 
tunneling experiments have the characteristics of an 
electron-bosonic mode coupling. The "41 mev" spin res- 
onance mode, a prominent feature in the spin excitation 
spectrum, is a natural candidate for electrons to cou- 
ple to |^,^|^. However, there are also phonon modes of 
similar energies, and they may instead have the strongest 
influence on the single electron spectra 

HUSH At 

present, ARPES alone appears inadequate to differen- 
tiate the two scenarios. Sometime ago, several of us 
proposed [13 an FT-IETS STM technique as a compli- 
mentary means to study this issue. The technique takes 
advantage of the pioneering work of the Fourier trans- 
formed STM IllL |l3| , and combines it with the vintage 
lETS [UlIlEl- 

Central to this technique is the Fourier transform of the 
energy derivative of tunneling conductance map in real 
space (Pl/dV'^{v,eV) I / dV^ {q, eV) . This q-space 
map, which can also be called Fourier map, contains in- 
formation about inelastic scattering in the system. Theo- 
retically, one finds peaks in q space and energy eV in this 
Fourier map of lETS that are related to the inelastic scat- 
tering off some collective excitations in the system. In the 
case of electron-spin mode coupling, the FT-LDOS near 
an ordinary potential scattering center (a impurity in 
the Nambu space) at the energy oi E = ±(Ao -|- flo) was 
shown to have sharp features at momenta close to (tt, tt). 
(Here Ap is the maximum of the d-wave superconducting 
energy gap and fig the mode energy.) 

Recently preliminary results from the first FT-IETS 



STM experiment has been reported in BSCCO [lij . 
While features are observed in the Fourier-transformed 
d^I/dV^ at the expected energy range, observed inten- 
sity near (tt, tt) is low. Instead, the strongest intensity 
appears at the wavevectors parallel to the Cu-0 bond 
directions. These experimental results have in turn mo- 
tivated us to compare the FT-IETS spectra near a po- 
tential scatterer in the cases of electrons coupled to the 
spin resonance and various phonon modes |l7l |. It was 
shown that all cases contain sharp features near (tTjTt), 
in disagreement with the experimental spectrum. 

This raises the question of whether the scatterer 
correctly describes the impurities in Bi2Sr2CaCu208+5 
(BSCCO). The origin of ubiquitously observed nanoscale 
inhomogeneity in BSCCO has been investigated by 
studying the correlation between the inhomogeneity and 
the position of oxygen dopants j^^. It is shown that 
the local electronic states are not associated with charge 
density variations. To account for those features, Nun- 
ner et al. proposed to look at the effect of random 
pairing potential fluctuations, the so-called n disorder. 
In what follows, we will address the effect of ti disorder 
on FT-IETS signatures. 

In this Letter, we take the disorder in pair potential 
as the scattering center, which is of ri character in the 
Nambu space, and study the Fourier component of the 
energy derivative local density of states, d'^I/dV'^, at 
E « ±(rio + Aq) around such a scatterer. We consid- 
ered a few typical bosonic modes as a possible scatter- 
ing modes that produce lETS fingerprints. We find that 
the results for ri are qualitatively different from the case 
of potential disorder: 1) there are no strong signatures 
near in the q-space near (tt, tt) in any of the electron- 
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boson couplings; 2) the highly anisotropic coupling of 
electrons to the out-of-plane out-of-phase oxygen buck- 
ling Big phonon mode, gives rise to a Fourier pattern 
similar to the lETS-STM experiment in BSCCO 
Our results are also consistent with the in-plane breath- 
ing mode although the agreement with the data is not as 
good. 

We start with a model Hamiltonian for a two- 
dimensional d-wave superconductor with the coupling of 
electrons to bosonic modes: 



i-th site, summed over two spin components: 
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Here the bare BCS Hamiltonian, 
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(1) 

-kiCkt), where 



the normal state energy dispersion is given by [2C 



4io cos k^. cos fc,. 



= —2ti(cOS k-j: -t- cos ky) — '±L2 <^^o ^^•=> '"y 

2i3(cos2fca; + cos2ky) — 4<4(cos2fc^cosfcy -t- 
cos fca; cos 2fcy) — 4t5 cos 2fc2, cos 2fcy — n, with ti ~ 1, 
t2 = -0.2749, t3 = 0.0872, ti = 0.0938, ts = -0.0857, 
and yU = —0.8772, and the d-wave gap disper- 
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sion Ak = 
fied explicitly, 
of ti hereafter, 
to bosonic modes 
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speci- 

encrgy is measured in units 
The coupling of the electrons 
is modeled by the Hamiltonian 

.9>^(k,q)4+q,^Cka(&<.q + 6t,-q) 

Big {v = 1) and the in- 
breathing {v ~ 2) modes, while 



for 



2JVr '^k+q,(T(^q ■ ''■(T(T')Ck,, 

the spin resonance mode. For the phonon modes, we 
consider the cases where the coupling matrix clement is 
either highly anisotropic, dependent on both k and q, or 
is only q dependent. In the following we use the notation 
Big-1 and br-1 for the former type of phonon modes 
while Big-11 and hr-l\ for the latter type. Detailed 
expression of the coupling matrix elements for these 
types of phonon modes can be found in Ref. The 
third term describes the quasiparticles scattered off a 
Ti impurity due to the inhomogeneity in pair potential 
rather than off a conventional T3 impurity. In the 
following, we consider a single ti impurity in a d — wave 
superconductor. The resulting Fourier pattern should 
survive a white-noise random distribution of such ti 
impurities in a realistic system. The impurity part of 
Hamiltonian can then be written as: 
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H.C.] , (2) 



where r]s = 1(^1) for S = S:(y). 

To be relevant to recent experimental realization, 
where no impurity-induced resonance state was observed, 
we assume the n impurity to have a weak scattering po- 
tential SA. In this limit, we employ the Born approxi- 
mation and arrive at the correction to the LDOS at the 



Sp{n,E) = Vry5lm[g(i,0;£- 

TT ^ — ' 



'ij)fiQ{S,i;E+i-f)]ii 



(3) 

where Q is the Green's function dressed with the 
bosonic renormalization effect and defined in the Nambu 
space 01 . From the perspective of the lETS, the energy 
derivative of the LDOS, is the quantity we are 

most interested in. It corresponds to the derivative of the 
local differential tunneling conductance (i.e., (Pl/dV^). 
The Fourier component of the differential LDOS is then 



dE 



e with the spectral 



d(5p(q,_E) 
dE 



given by ' ''^^ ^ E 
weight defined as P(q, E) 

We consider here for comparison the coupling of elec- 
trons to spin resonance mode. Big and breathing phonon 
modes. For the numerical calculation, we take the super- 
conducting energy gap Aq = 0.1, the frequency of all 
bosonic modes VLq = 0.15. The ti impurity scattering 
strength 6 A is taken to be 50% of the superconducting 
energy gap. The coupling strength for all types of bosonic 
modes is calibrated to give at the Fermi energy E = Q 
an identical frequency renormalization factor in the self 
energy. The same procedure as in Ref. 0| is followed to 
obtain the Fourier spectral weight P(q, E). 

In Figs. ^ we present the results of the Fourier spec- 
trum, P{q,E), at the energy E = — (Ao + Hq) for a 
d-wave superconductor with the electronic coupling to 
the bosonic modes. For comparison, the same spectrum 
is also shown (last panel) for the case of no mode cou- 
pling. Note that the case without the mode couphng, the 
energy Qq has no special meaning in the context of the 
electronic properties, and the energy E ~ — (Aq + flo) is 
chosen merely for comparison to the case of mode cou- 
pling. The energy E — — (Aq + ^lo) corresponds to the 
position where the bosonic modes are excited, signalin 
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a peak in the lETS d I/dV -V tunneling spectrum 
First of all, the Fourier maps for all cases does not dis- 
play any peak structure at large q near (±7r, ±7r) and 
(±7r, TJlLl which appears persistently with the T3 scat- 
tering [l2|. Instead the Fourier spectral weight is min- 
imal in intensity (dark blue area in the figure) in these 
regions. The map for the case of electronic coupling to 
the Big-l phonon mode shows locally strongest intensity 
(red spots) at q about (±^,0) and (0, ±^). The peak 
intensity at q near (±^, ±^) and (±^, T^) much 
weaker than those along the bond directions. The map 
for the case of the coupling to the br-l phonon mode ex- 
hibits locally strongest intensity (red spots) at q near 
(±^,0) and (0, ±1^). In addition, each of these red 
spot has a double-tail structure, which is absent in the 
case of Big-l mode coupling. The maps for the cases 
of the Big-ll and spin resonance mode coupling exhibit 
similar q structure. The finite intensity is uniformly dis- 
tributed on a circular strip near |q| = ^ and becomes 
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FIG. 1: The Fourier spectrum of the energy derivative LDOS 
at the energy E = — (Ao + f2o) for a d-wave superconductor 
with the electronic coupling to the Big-1, br-1, Big-11, br-ll, 
spin resonace modes. For comparison, the spectrum is also 
shown for the case of no mode coupling. 



stronger as q approaches the zero point. No locally 
distinguishable strongest intensity peak can be identi- 
fied at q = (±^,0) and (0,±3f). The map for the 
coupling to the br-ll phonon mode exhibits locally the 
highest intensity (red spots) at q near (±^,±^) and 
(±^, =F^)- No peaks are found at q near (±^,0) and 
(0, ±^). The map for the case of no mode coupling 
shows an eight-tail star shape at q = (0,0), which con- 
sists of the head-on overlap of four red spots, such as 
those appearing in the case of the br-1 phonon mode 
coupling each with two tails. As we have already em- 
phasized, experimentally, the Fourier map of (Pl/dV^ 
shows intensity peaks only at q = (±^,0) ± 15% and 
(0,±^) ± 15% 113. Therefore, by comparison with the 
experimental data, our new FT-IETS STM analysis also 
supports the notion that the elec- 

tronic band must be renormalized by its coupling to the 
bosonic modes. In particular, the results based on the 
scenario of highly anisotropic coupling of electrons to the 







FIG. 2: The Fourier spectrum of the derivative of the LDOS 
is shown at the various values of the energy for the case of 
the electronic coupling to the -B19-I phonon mode. Here the 
energy has been measured by scaling Ao — 30 meV. 



In Fig. 12 we present the energy evolution Fourier pat- 
tern for the electronic coupling to the Big-1 mode. It 
shows that the intensity peak structure along the bond 
direction is robust against the energy change. The char- 
acteristic q vector, at which the locally highest intensity 
is located, decreases slightly with the increased energy. 
This result is also not inconsistent with the preliminary 
experiment. 

Our results demonstrate the important role the charac- 
ter of the scattering center plays in the Fourier spectrum. 
To explore this issue further, we note that, in the case of 
Ti scattering considered here, the Fourier spectrum can 
be expressed as follows. 



<^p(q; E)^ ^ Ek(cos K - cos ky) 

x{[A{k; E)K{k + q;E)+ A'(k; E)A{k + q; E)] 

+ [J(k; E)B{k + q;E) + B{k; E)J{k + q; E)]} , (4) 



where 



A{k;E) lm[gn{k;E + tj)] , (5) 

B{k;E) = Rc[gu{k;E + tj)] , (6) 

J{k;E) = --lm[gi2{'k;E + ij)], (7) 

TT 

K{k;E) = Rc[g,2{k;E + tj)] . (8) 



This expression shows that, for the ti scattering, the 
Fourier spectrum is determined by the k-summation of 
the product terms constituting the imaginary (real) parts 
of the single-particle (^/n) with the real (imaginary) 
parts of anomalous (t/i2) Green's function in the super- 
conducting state, weighted by a d-wave type form fac- 
Big phonon mode arc in best agreement with the lETS- tor coskx — cos ky. This scattering process with the 
STM data in BSCCO. ti impurity is significantly different than the case of a 
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Ta impurity scattering |l7j . where the convolution takes 
place between the real and imaginary parts of the same 
Green's function without the form factor. This differ- 
ence of the scattering process matters significantly in the 
resulting Fourier map. As shown in Eq. |0J, the form fac- 
tor cos kx — cos ky appearing in the ri scattering case is 
identically zero along the diagonals in the first Brillouin, 
but reaches a maximum at the AI points [k = (±7r, 0) 
and (0, ±7r)] on the zone boundary . It then follows that 
any stronger intensity from the product of AK , BJ con- 
nected by a q oriented close to the diagonals is strongly 
suppressed, while the intensity connected by q oriented 
parallel to the bond direction is enhanced. For the elec- 
tronic coupling to the Big-I phonon mode, it has been 
found that there are moderately strong intensity on 
the two split beams oriented perpendicular to the zone 
boundary at M points in the function A, B, J, and K. 
The form factor cos k^ — cos ky then tips the relative con- 
tribution from the product AK and BJ , giving rise to 



locally highest intensity at q ■ 



,0) and (0,±^) in 



the Fourier map (see the first panel of Fig. P). These split 
beams of intensity are absent for the electronic coupling 
to other modes. 

Our results naturally suggest additional means to fur- 
ther explore the electron-bosonic mode coupling exper- 
imentally. For instance, a Zn impurity acts as a non- 
mangetic potential center - & scatterer. In this case, 
the sharp features near (tt, tt) should be observed in the 
FT-IETS spectrum. In the case of low-energy elastic 
scattering interference of quasiparticles, related effects 
have in fact been demonstrated. Indeed, strong signa- 
tures near (±7r, ±7r) and (±7r, =F7r) appear in the theo- 
retical spectra near a T3 scatterer |2l|. These features 
arc not observed experimentally in the stoichiomctrical 
BSCCO [ll[l3, but arc seen around a nomagnctic Zn 
impurity in the doped BSCCO. 

In conclusion, we have studied, for the first time, 
the Ti-impurity induced Fourier pattern of the energy 
derivative local density of states in a d-wave supercon- 
ductor with the coupling of electrons to various bosonic 
modes. We consider -Big, half-breathing, and spin (tt, tt) 
modes. Our results show that, at the mode energy 
shifted by the d-wave superconducting gap energy Ag, 
i.e., E = ±(Ao -I- flo), the coupHng of electrons to the 
Big or breathing phonon modes, gives rise to a Fourier 
pattern similar to the preliminary Fourier transformed 
lETS-STM experiment in BSCCO The coupUng of 
electrons to the spin resonance mode, on the other hand, 
yields a Fourier spectrum that is inconsistent with the 
experiment. These results do not necessarily rule out 
the role of the spin-spin interactions as being relevant for 
superconductivity in BSCCO, instead they imply that 
electron-phonon coupling has a strong impact on the su- 
perconducting electronic structure. These results have 
important implications for our understanding of the elec- 
tronic properties of the cupratcs. They also demonstrate 



the potential of the FT-IETS STM technique, and high- 
light the imp ortance of ri scattering in the impurity-free 

BSCCO Ha. 
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